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Abstract 

We extend the application of the techniques developed within the framework of the 
pseudo-Hermitian quantum mechanics to study a unitary quantum system described by 
an imaginary PT-symmetric potential v{x) having a continuous real spectrum. For this 
potential that has recently been used, in the context of optical potentials, for modelling the 
propagation of electromagnetic waves travelling in a wave guide half and half filed with 
gain and absorbing media, we give a perturbative construction of the physical Hilbert 
space, observables, localized states, and the equivalent Hermitian Hamiltonian. Ignoring 
terms of order three or higher in the non-Hermiticity parameter £, we show that the 
equivalent Hermitian Hamiltonian has the form ^— + ^- Y^=®{ a n{ x )-,P 2n } with a n (x) 
vanishing outside an interval that is three times larger than the support of v(x), i.e., in 
2/3 of the physical interaction region the potential v(x) vanishes identically. We provide 
a physical interpretation for this unusual behavior and comment on the classical limit of 
the system. 

PACS numbers: 03.65.-w 
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1 Introduction 



During the past seven years there have appeared over two hundred research papers on VT- 
symmetric quantum systems. This was initially triggered by the surprising observation of Bessis 
and Zinn- Justin and its subsequent numerical verification by Bender and his co-workers j]Q that 
certain non-Hermitian but PT-symmetric Hamiltonians, such as 

H = p 2 + x 2 + iex 3 with e 6 IR + , (1) 

have a purely real spectrum. This observation suggested the possibility to use these Hamiltoni- 
ans in the description of certain quantum systems. Since the PT-symmetry of a non-Hermitian 
Hamiltonian H, i.e., the condition [H,VT] = 0, did not ensure the reality of its spectrum, a 
crucial task was to seek for the necessary and sufficient conditions for the reality of the spec- 
trum of a given non-Hermitian Hamiltonian H . This was achieved in [2] where it was shown, 
under the assumptions of the diagonalizability of H and discreteness of its spectrum, that 
the reality of the spectrum was equivalent to the existence of a positive-definite inner prod- 
uct (•,•)+ that rendered the Hamiltonian self-adjoint, i.e., for any pair of state vectors 
(1>,H<p) + = (Hil>, <!>)+. 

Another condition that is equivalent to the reality of the spectrum of H is that it can 
be mapped to a Hermitian Hamiltonian h via a similarity transformation [21 E|; there is an 
invertible Hermitian operator p such that 

H = p~ 1 hp. (2) 

The positive-definite inner product (■,■)+ and the operator p entering (J2J) are determined by a 
positive-definite operator i] + according to [21 E] 

(-,-)+ :=<-kf>, (3) 
P = VvT, (4) 

and the Hamiltonian satisfies the r/ + -pseudo-Hermiticity condition j3J: 

= V+HV+ 1 - (5) 

Here (-|-) stands for the standard (L 2 ) inner product that determines the (reference) Hilbert 
space 7i as well as the adjoint W of H, J5]. 1 

It is this, so-called metric operator, i] + that determines the kinematic structure (the physical 
Hilbert space and the observables) of the desired quantum system. Note however that r] + is not 
unique [HI El IE]- 2 In [2] we have not only established the existence of a positive definite metric 
operator r/ + and the corresponding positive-definite inner product (•,•)+ for a diagonalizable 
Hamiltonian 3 with a discrete real spectrum, but we have also explained the role of antilinear 

1 The adjoint At of an operator A is the unique operator satisfying, for all ^i^gH, (ip\A'(f>) = (Atp\(f>). A 
is called Hermitian if A^ = A. 

2 it is only unique up to symmetries of the Hamiltonian, [7]. 

3 For a treatment of non-diagonalizable pseudo-Hermitian Hamiltonians see [H1E1E]- Note that diagonal- 
izability of the Hamiltonian is a necessary condition for applicability of the standard quantum measurement 
theory . It is also necessary for the unitarity of the time-evolution, for a non-diagonalizable Hamiltonian is 
never Hermitian (its evolution operator is never unitary with respect to a positive-definite inner product, 

dnni. 
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symmetries such as VT and offered a method for computing the most general r] + . An alternative 
approach that yields a positive-definite inner product for a class of PT-symmetric models is 
that of [12 . As shown in the CVT- inner product proposed in ^2] is identical with the 

inner product (■, ■)+ = (-|?y+-) for a particular choice of 77+. 

Under the above mentioned conditions every Hamiltonian having a real spectrum determines 
a set Uh+ of positive-definite metric operators. To formulate a consistent unitary quantum 
theory having H as its Hamiltonian, one needs to choose an element r] + of Uh+- a Each choice 
fixes a positive-definite inner product (•,■)+ and defines the physical Hilbert space 7i p h y s and 
the observables. The latter are by definition jT^j the operators O that are self-adjoint with 
respect to (•,•) + , alternatively they are ?7 + -pseudo-Hermitian. These can be constructed from 
Hermitian operators o acting in 7i according to [5] 

= p- 1 op. (6) 

In particular, one can define 77+-pseudo-Hermitian position X and momentum P operators 
[Tol lo]. express if as a function of X and P, and determine the underlying classical Hamiltonian 
for the system by letting h — > in the latter expression, Alternatively, one may calculate 

the equivalent Hermitian Hamiltonian h and obtain its classical limit (again by letting h — > 0). 

Another application of the 7? + -pseudo- Hermitian position operator X is in the construction 
of the physical localized states: 

\Z {X) ) := p-'lx). (7) 

These in turn define the physical position wave function, ty(x) := (^ x \ip) + = (x\p\^), and the 
invariant probability density, 

oh) - mx)l * = »^»^» 2 (8) 

for a given state vector \ip), [5| ITK]. 

The above prescription for treating PT-symmetric and more generally pseudo- Hermitian 
Hamiltonians with a real spectrum has been successfully applied in the study of the VT- 
symmetric square well in and the cubic anharmonic oscillator (P) in |T6J. 5 Both these 
systems have a discrete nondegenerate energy spectrum, and the results of jH |2] are known to 
apply to them. The aim of the present paper is to seek whether these results (in particular 
the construction method for r] + ) may be used for treating a system with a continuous spec- 
trum. 6 This question is motivated by the desire to understand field-theoretical analogues of 
■pT-symmetric systems which should admit an S'-matrix formulation. Furthermore, there are 
some basic questions related to the nonlocal nature of the Hermitian Hamiltonian h and the 
pseudo-Hermitian observables such as X and P especially for PT-symmetric potentials with a 
compact support (i.e., potentials vanishing outside a compact region). 

To achieve this aim we will focus our attention on a simple toy model recently considered as 
an effective model arising in the treatment of the electromagnetic waves travelling in a planar 

Alternatively one may choose sufficiently many operators with real spectrum to construct a so-called irre- 
ducible set of observables which subsequently fixes a metric operator r)+, |14j . 
5 See also Q7]. 

6 The question whether the theory of pseudo-Hermitian operators as outlined in 0J|2] is capable of treating a 
system having scattering states was posed to the author by Zafar Ahmed during the 2nd International Workshop 
on Pseudo-Hermitian Hamiltonians in Quantum Physics, held in Prague, June 14-16, 2004. 
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slab waveguide that is half and half filed with gain and absorbing media, This model has 
a standard Hamiltonian, 

H = t- + „(*), (9) 
and a "PT-symmetric imaginary potential, 

{0 for \x\ > ^ or x = 
< for ie(-|,0) (10) 
-iC for x G (0,|), 

where L G (0, oo) is a length scale, £ G [0, oo) determines the degree of non-Hermiticity of the 
system, and 9 is the step function: 

{0 for x < 
\ for x = (11) 
1 for x > 0. 

The Hamiltonian (jHJ) differs from a free particle Hamiltonian only within (—■§,■§) where it 
coincides with the Hamiltonian for the "PT-symmetric square well [TTH Ej . 

It is important to note that unlike in JTHj we will consider the potential (jlOj) as defining a 
fundamental (non-effective) quantum system having a unitary time-evolution (and S- matrix). 
Therefore our approach will be completely different from that pursued in [IB] and the earlier 
studies of effective (optical) non-Hermitian Hamiltonians, [20]. 

To the best of author's knowledge, the only other non-Hermitian Hamiltonian with a con- 
tinuous (and doubly degenerate) spectrum that is shown to admit a similar treatment is the 
one arising in the two-component formulation of the free Klein-Gordon equation • Com- 

pared to this Hamiltonian defines a technically much simpler system to handle, because 
it is essentially a tensor product of an ordinary Hermitian Hamiltonian and a 2 x 2 matrix 
pseudo-Hermitian Hamiltonian. 



2 Metric Operator 

The essential ingredient of our approach is the metric operator rj + . For a diagonalizable Hamil- 
tonian with a discrete spectrum it can be expressed as 

da 

i]+ = X^SZ l</>n,a)(0 n ,a|, (12) 

n a=l 

where n, a, and d n are a spectral label, a degeneracy label, and the multiplicity (degree of 
degeneracy) for the eigenvalue E n of H, respectively, and {\(p n ,a)} is a complete set of eigen- 
vectors of that together with the eigenvectors \ip n ,a) of H form a biorthonormal system, 

Now, consider a diagonalizable Hamiltonian with a purely continuous doubly degenerate real 
spectrum {Ek}, where k G (0, oo). We will extend the application of (|12|) to this Hamiltonian 
by changing • • • to f dk - • ■ . This yields 

POD 

T] + = / dfc(|0fc,+)<0fc,+| + |0fc,-><0 fc ,-|), (13) 
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where we have used ± as the values of the degeneracy label a, [22] • The biorthonormal system 
{\ip k ,a),\4> k ,a)} satisfies 

H\ip k ,a} = E k \ip k ,a), H j \<p k , a) = E k \<fr k , a), (14) 

POO 

{(j) k ,a\ip e ,b) = 5 ab 5(k-£), / (\if) k ,+)(<f> h ,+\ + \ij; k ,-)(<p k ,-\)dk=l, (15) 



where 5 a b and S(k) stand for the Kronecker and Dirac delta functions, respectively, k G (0, oo), 
and a, b G {— , +}, 

We define the eigenvalue problem for the Hamiltonian Q using the oscillating (plane wave) 
boundary conditions at x — ±oo similarly to the free particle case which corresponds to ( = 0. 
To simplify the calculation of the eigenvectors we first introduce the following dimensionless 
quantities. 

x :=(£)*, p:=(A)p, Z:=(#R, H := H = p 2 + v(x), (16) 

{0 for |x| > 1 or x = 
iZ for xG(-l,0) (17) 
-iZ for xG(0, 1). 

The eigenvalue problem for the scaled Hamiltonian H corresponds to the solution of the 
differential equation 

d 2 

+ vfx) - E fc 



dx 2 



ifj(x) = 0, (18) 



that is subject to the condition that if) is a different iable function at the discontinuities x = 
— 1,0,1 of v. Introducing if>x : (— oo, — 1] — > C, if)- : [—1,0] — > C, if> + : [0,1] — > C, and 
ip2 '■ [1, oo) — > C according to 



i\)\ (x) for x G (— oo, — 1] 

* ip+(x) for x G [0, 1] 
tp2 (x) for x G [1, oo) 



we have 



^i(-i) = V-(-i), V>l(-i) = ^(-i), (20) 

^-(0) = v+(o), ^_(0)=^(0), (21) 

V + (l) = ^2(1), V;(l)=^(l)- (22) 

Now, imposing the plane-wave boundary condition at x = ±oo and demanding that the eigen- 
functions if; be "PT-invariant, which implies 

^_(0)=^ + (0)*, V-(0) = -V4(0), (23) 

we find E k = k 2 , i.e., the spectrum is real positive and continuous, and 

^i(x) = A ie ikx + B ie ~ ikx , if) 2 (*) = A 2 e ikx + B 2 e~ ikx , if)±{x) = A±e ik±x + B ± e' ik±x , (24) 
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where 



k± := Vk 2 ± iZ, 

ik 

A, — At — —^= 



ik 



Aj 



2vr 
cos 



[L_(ife)u + v] 
ik- sin 



57T 



k 
±1/2 




(25) 

+ (26) 



i sin 



(27) 
(28) 



and u,i;Gl are arbitrary constants (possibly depending on k and/or Z and not both vanish- 
ing). 

The presence of the free parameters u and v is an indication of a double degeneracy of 
the eigenvalues = k 2 . We will select u and v in such a way as to ensure that in the limit 
Z — > we recover the plane- wave solutions of the free particle Hamiltonian, i.e., we demand 



(29) 



lim.z_>o "0 ( x ) = e /y2n. This condition is satisfied if we set 



u 



±1. 



In the following we use the superscript ± to identify the value of a quantity obtained by setting 
u — 1 and v = ±1. In this way we introduce A 1 , Bf, Af, B^-, A±, B±, and ip ± . The latter 
define the basis (generalized eigenvectors \ipk, ±) by (x|0fc, ±) := V ;± ( x )- 

The next step is to obtain |0^,±). In view of the identity IT" = H\z-*-z, we can easily 
obtain the expression for the eigenfunctions (ft of W . Introducing 



x 



4>i (x) for x G (-oo,-l] 

0_(x) for xG[-l,0] 

0+(x) for x G [0, 1] 

2 (x) for x G [1, oo), 



(30) 



we have 

0i (x) = C ie ifex + Die- te , 
where 



C 2 e 



£> 2 e 



-ikx 



Ci 



/27T 

L+(*) :=L_(-*)* 



2tt 



[L + (-k)r + K + (-k)s] 



in 



r + 



=Fl/2 



:= -K_(-kY 



D-t 



in 



=Fl/2 



(31) 

(32) 
(33) 

(34) 



and r, s G M are (possibly k- and/or Z-dependent) parameters that are to be fixed by imposing 
the biorthonormality condition ()15)) . The latter is equivalent to a set of four (complex) equations 
(corresponding to the four possible choice for the pair of indices (a, b) in the first equation in 
(115(0 which are to be solved for the two real unknowns r and s. This together with the presence 
of the delta function in two of these equation make the existence of a solution quite nontrivial. 
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We checked these equations by expanding all the quantities in powers of the non-Hermiticity 
parameter Z up to (but not including) terms of order two and found after a long and tedious 
calculation (partly done using Mathematica) that indeed all four of these equations are satisfied, 
if we set r = u = 1 and s — v — ±1. Again we will refer to this choice using superscript ±. In 
particular, we have ^ = ip^z^-z and (x\(pk, ±) := ^(x). 

Having obtained |0fc,±) we are in a position to calculate the metric operator (fT3"j) . We 
carried out this calculation using first order perturbation theory in Z. It involved expanding 
the 4>f (x) , ip^ ( x ) ? an d ^±( x ) m powers of Z, substituting the result in 

POO 

<x|77 + |y> = / [0+(y)* + (x) + </>-(yT <p-(*)]dk (35) 
Jo 

which follows from (|13)1. and using the identities: 

/oo /* oo iak f oo %ak ibk 

e iak dk = 2n5(a), / ^— dk = in sign(a), / - — dA; = tt(|6| - |a|) (36) 
oo J —oo ™ J —oo 

(where a, 6 G K and sign(a) := 0(a) — 0(— a)) to perform the integral over /c for all 16 possibilities 
for the range of values of the pair of independent variables (x, y) in (j35[). To simplify the 



presentation of the result, we introduce: I\ := (— oo, — 1), I : = (—1,0), 7 + := (0,1), I 2 := 
(1, 00) and define the functions E M) „ : 1^ x /„ — > C by 

^, v (x,y) := (x-\v+\y) for all xel^jEh, fi, v G {1, -, +, 2}. 
Then after a very long calculation we find 

£i,l(x,y) = o(x-y) + §sign(x-y)Z + 0(Z 2 ), (37) 

£_,i(x,y) = |(2-x-y-|x + y + 2|)Z + 0(Z 2 ), (38) 

f+,i(x,y) = i(2-x-y-|x + y + 2|)Z + 0(Z 2 ), (39) 

£ 2)1 (x,y) = |(4 + 2|x + y|-|x + y + 2|-|x + y-2|)Z + 0(Z 2 ), (40) 

£i,_(x,y) = -|(2-x-y-|x + y + 2|)Z + 0(Z 2 ), (41) 

f_,_(x,y) = 0(x-y)-|sign(x-y)(x + y)Z + O(Z 2 ), (42) 

£ + ,-(x,y) = ||x + y|Z + 0(Z 2 ), (43) 

£ 2 ,_(x,y) = |(2 + x + y-|x + y-2|)Z + 0(Z 2 ), (44) 

£i )+ (x,y) = -|(2-x-y-|x + y + 2|)Z + 0(Z 2 ), (45) 

f-.+(x,y) = -||x + y|Z + 0(Z 2 ), (46) 

£+,+(x,y) = 0(x-y) + isign(x-y)(x + y)Z + O(Z 2 ), (47) 

£ 2)+ (x,y) = |(2 + x + y-|x + y-2|)Z + 0(Z 2 ), (48) 

£ 1)2 (x,y) = -|(4 + 2|x + y|-|x + y + 2|-|x + y-2|)Z + 0(Z 2 ), (49) 

£_, 2 (x,y) = -f(2 + x + y-|x + y-2|)Z + 0(Z 2 ), (50) 

£+, 2 (x,y) = -|(2 + x + y-|x + y-2|)Z + 0(Z 2 ), (51) 

£ 2 , 2 (x,y) = 0(x-y) + §sign(x-y)Z + O(Z 2 ), (52) 



where 0(Z 2 ) stands for terms of order two and higher in powers of Z. It is quite remarkable 
that we can obtain from ()37j) - f!52l) a single formula for (x|r/ + |y) which is valid for all x, y G R, 
namely 

<x|77+|y> = 5(x - y) + § (4 + 2|x + y| - |x + y + 2| - |x + y - 2|) sign(x - y) Z + G(Z 2 ). (53) 
Note that (x|?y + |y)* = (j\i] + \x) which is consistent with the Hermiticity of r) + . 
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3 Physical Observables and Localized States 

The physical observables of the system described by the Hamiltonian @ are obtained from the 
Hermitian operators acting in 7i = L 2 (M.) by the similarity transformation ©. This equation 
involves the positive square root p of r] + which takes the form jTH| 

p ±l = e^l\ (54) 

if we express i] in the exponential form 

r] + = e~ Q . (55) 
In view of (154)1 and the Backer-Campbell-Hausdorff identity, 

e~ A B e A = B + [B, A] + ± [[B, A],A]-\ (56) 

(where A and B are linear operators), physical observables (jOJ) satisfy [TH] : 

= o-^[o,Q) + ^[[o,Q),Q}+--- . (57) 

If we expand r/ + and Q in powers of Z, 

oo oo 

V+ = l + J2v +e Z e , Q = Y,QeZ e , (58) 
i=i e=i 

where rj +i and Qe are Z-independent Hermitian operators, we find using that 

Qi = -?7+i> = -v+2 + \ r i+i- ( 59 ) 

Combining this relation with ([57)1 . we have 

O = o - i[o, Qi] Z + i(-4[o, Q 2 ] + [[o, gj, Ox]) Z 2 + 0(Z 3 ). (60) 

In the following we calculate the r/ + -pseudo-Hermitian position (X) and momentum (P) 
operators, JH]; up to (but no including) terms of order Z 2 . This is because so far we have only 
calculated r/ +1 which in view of flB^j) satisfies 

(x|77 +1 |y> = |(4 + 2|x + y| - |x + y + 2| - |x + y - 2|) sign(x - y), VxjGR. (61) 

Substituting the scaled position (x) and momentum (p) operator for o in (|60|) . using (jHTjl . and 
doing the necessary algebra, we find 

(x|X|y) = x5(x-y) + ^(4 + 2|x + y|-|x + y + 2|-|x + y-2|)|x-y|Z + 0(Z 2 ), (62) 
( x l p |y) = ~id x S(x - y) + | [2sign(x + y) - sign(x + y + 2) 

-sign(x + y - 2)] sign(x -y)Z + 0(Z 2 ), (63) 

where X := 2X/L and P := LP/(2h) are dimensionless ^ + -pseudo-Hermitian position and 
momentum operators, respectively. 
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As seen from f!62|) . both X and P are manifestly nonlocal and non-Hermitian (but pseudo- 
Hermitian) operators. Furthermore, 

(x|P|y) = (x|p|y> + 0{Z 2 ) for x $ [-1, 1] and y $ [-1, 1]. 

If we scale back the relevant quantities in ()62|) and (jp^|) according to (fH)j) . we find 

(x\X\y) = x 5{x-y) + i${2L + 2\x + y\-\x + y + L\-\x + y-L\)\x-y\( + 0{( 2 ),{te) 
(x\P\y) = -ihd x 6(x - y) + f| [2 sign(x + y) - sign(a; + ?/ + L) 

-sign(x + y-L)] sign(x - y) ( + 0(C)- (65) 

Again note that the contributions of order £ to P vanish, if both x and y take values outside 

L 2 ' 2 J • 

Next, we compute the localized states t^ x > of the system. The corresponding state vectors 
are defined by (J7J). Using this equation as well as (p>4"j) . (p)%j) . (joTlj) . (jBTj) . and (fTo^l we have the 
following expression for the x-representation of a localized state £^ centered at y G KL 

<x|£<*>> = 5(x - y) - ^ (2L + 2\x + y\ - \x + y + L\ - \x + y - 1 1) sign(x -y) + <D(( 2 ). (66) 

Because the linear term in ( is imaginary, the presence of a weak non-Hermiticity only modifies 
the usual (Hermitian) localized states by making them complex (non-real) while keeping their 
real part intact. Note however that for a fixed y the imaginary part of does not 

tend to zero as \x — y\ — > oo. This observation which seems to be in conflict with the usual 
notion of localizability has a simple explanation. Because the usual x operator is no longer 
an observable, it does not describe the position of the particle. This is done by the pseudo- 
Hermitian position operator X; it is the physical position wave function := (£ {x) ,^)+ that 
defines the probability density of localization in space (jHJ). The physical position wave function 
for the localized state £W is given by (^ x \ = (x\y) = 5(x — y) which is the expected 

result. 

In summary, the notion of localizability in space is directly linked with the choice of the 
physical position operator. An important by-product of the recent intensive investigation of 
non-Hermitian PT-symmetric systems is the realization of the fact that one my formulate a 
unitary quantum system for which the choice of the Hilbert space and observables, particularly 
the position operator, is not a priori fixed. 



4 Equivalent Hermitian Hamiltonian and Classical Limit 

The calculation of the equivalent Hermitian Hamiltonian h for the Hamiltonian Q is similar 
to that of the physical observables. In view of (J2J), (|54|). (|5b|). (|5S]). and the last equation in 
(116)) which we express as 

H = p 2 + iz/(x)Z with i/(x) :=0(x + l) + 0(x-l)- 20(x), (67) 

we have 

h = p 2 + h 1 Z + h 2 Z 2 + 0(Z 3 ), (68) 

hi := zz/(x) + ^[p 2 ,Qi], (69) 

h 2 : = ^{4[p 2 ,g 2 ] + 4^Kx),gi] + [[ P 2 ,gi],gi]}. (70) 
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where 

h:= pRp- 1 = mL 2 h/(2h 2 ) (71) 

is the dimensionless Hermitian Hamiltonian associated with H. 
Next, we substitute and (JoTj) in the identity 

(x|[ P 2 ) g 1 ]|v) = (9 y 2 -9 x 2 )(x|g 1 |y), 

and perform the necessary algebra. We then find (x|[p 2 , Qi]|v) = — 2iu(x)5(x — y). Therefore, 

\p 2 ,Q 1 )=-2iu(x), (72) 

and in view of ()69|) 

hi = 0. (73) 

This was actually to be expected, for both the operators appearing on the right-hand side 
of ()69|) are ant i- Hermitian, while its left-hand side is Hermitian. The fact that an explicit 
calculation of the right-hand side of (JBTJj) yields the desired result, namely ()73|) . is an important 
check on the validity of our calculation of rj +1 . It may also be viewed as an indication of the 
consistency and general applicability of our method, that was initially formulated for systems 
with a discrete spectrum Hn] . 
According to (JTHJ), 

h = p 2 + h 2 Z 2 + C(Z 3 ). (74) 

Hence, in order to obtain a better understanding of the nature of the system described by the 
Hamiltonian H, we need to calculate h 2 . Equations (|70|) and (J5HJ) suggest that this calculation 
demands a complete knowledge of 7] +2 which in turn requires the calculation of (x|?7 +2 |y) for 
all 16 possibilities for the ranges of x and y. This is an extremely lengthy calculation in which 
one must deal with infinite integrals of the form e mk / k n dk with n = 2, 3, 4. 7 We will 
not include the result of this calculation here, not only because it is too lengthy but most 
importantly because, as we will show in the following, the knowledge of (x|^ +1 |y) turns out to 
be sufficient for the calculation of h 2 . To see this we first employ ()72|) to express h 2 in the form 

h 2 = ^ (2[p 2 ,Q 2 ]+^Kx),g!]). (75) 

Now, we recall that p 2 , Q 2 , v{x) and Qi are all Hermitian operators. Therefore [p 2 ,Q 2 ] and 
z[i/(x),Qi] are respectively anti- Hermitian and Hermitian. In view of (|75|) and the Hermiticity 
of h 2 , this implies that 

[P 2 ,Q 2 ]=0. (76) 

Hence, 

h a = l[i/(x),g 1 ] = l[»7 +1 ,i/(x)], (77) 

where we have also made use of the first equation in ([59)1 . We should also mention that the 
identities (I72j) and (|75|) can be directly obtained from the pseudo-Hermit icity condition (jSJ) by 
substituting (J55j) in (JSJ) and using (J55j) an d (|58|). 

We can easily use (J61|) and (J77j) to yield the expression for the integral kernel of h 2 , namely 

(x|h 2 |y) = ^(4 + 2|x + y|-|x + y + 2| -|x + y-2|)sign(x-y)[i/(x)-Ky)], Vx,yGM. (78) 
7 These may be easily regularized as is well known in typical field theory calculations. 



10 



As seen from this equation, 

(x|h 2 |y) = 0, if x£ [-1,1] and y£ [-1,1]. (79) 

We can express h 2 as a function of the x and p by performing a Fourier transformation on 
the y variable appearing in (J7SJ), i.e., computing 

/oo 
(x|h 2 |y)e ip ^y. (80) 
-oo 

This yields h 2 as a function of x and p, if we order the factors by placing x's to the left of 
p's. We can easily do this by expanding (x|h 2 |p) in powers of p. Denoting the x-dependent 
coefficients by u n , we then have 

oo 

h 2 = 5>„(x)p", (81) 

n=0 

where we have made the implicit assumption that (x|h 2 |p) is a real-analytic function of p. 

The Fourier transform of (x|h 2 |y) can be performed explicitly. 8 We have instead used 
Mathematica to calculate (x|h 2 |p) and found the coefficients u n for n < 5. It turns out that 
indeed (x|h 2 |p) does not have a singularity at p = 0, and that co> , 1^2,^4 are real and vanish 
outside (—3,3) while uji,uj 3 ,lj 5 are imaginary and proportional to 9(x) — 1/2 outside (—3,3). 
As we will explain momentarily these properties are necessary to ensure the Hermiticity of h. 

Figures 1, 2 and 3 show the plots of real part of u n for n = 0, 2, 4 and the imaginary part 
of u n for n — 1, 3, 5. As seen from these figures (the absolute value of) u n sharply decreases 
with n, which suggests that a truncation of (J81|) yields a a good approximation for the action 
of h 2 on the wave functions with bounded and sufficiently small x-derivatives. 

If we use (p|h 2 |x) = (x|h 2 |p)* to determine the form of h 2 and suppose that ^ 2n (x) are real 
and c<j 2n+ i(x) are imaginary for all n = 0, 1, 2, 3, • • • , we find 

00 00 

h 2 = J2v n "n(x)* = $> 2 V> 2n (x) - P 2 " + W +1 (X)]. 
n=0 n=0 

Adding both sides of this relation to those of (JHTj) and diving by two, we obtain 



^ 00 

n=0 



(x), p 2 "}, a„(x) := w 2n (x) + iuj' 2n+1 (x), (82) 



where {•, •} stands for the anticommutator, a prime denotes a derivative, and we have made use 
of the identity: [/(x), p m ] = {if'(x), p m_1 }. It is important to note that because c<j 2n (x) are real 
and a; 2n+ i(x) are imaginary, a n (x) are real. Moreover, outside (—3,3), u; 2 „(x), u' 2n+1 (x), and 
consequently a n vanish. Therefore, we can express h 2 in the manifestly Hermitian form (|82j) 
with all the x-dependent coefficient functions vanishing outside (—3,3). Figure 4 shows the 
plots of a n for n = 0,1,2. They are all even functions of x with an amplitude of variations that 
decreases rapidly as n increases. 

8 One way of doing this is to use the integral representations of the absolute value and sign function, as given 
in J2HJ), to perform the y- integrations in (|80J) and use the identities 

/oo iau l l-oo p iauj ■ 

= ™ {e iak x) sign(a); / 6 au = g[i + {iak _ i )e **] sign(a); Va , k e M, 

to evaluate the remaining two integrals. The resulting expression is too lengthy and complicated to be presented 
here. 
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Figure 1: Graph of the real part of co (dashed curve) and u 2 (full curve). 
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Figure 2: Graph of the imaginary part of oj\ (dashed curve) and U3 (full curve). 
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Figure 3: Graph of the real part of uj^ (dashed curve) and the imaginary part 
of UJ5 (full curve). 



Next, we scale back the relevant quantities and use (fTfiJl . (|7TJ). (J74*)l . and (jH2| to obtain 

2 /-2 00 



a n (ar) := 2m (^) 



(83) 



n=0 



In view of the fact that a n and a n are real-valued even functions, h is a manifestly Hermitian 
V- and T-symmetric Hamiltonian. We can also express it in the form 



1 C 2 00 



34) 



n=2 



where 



7?) 



m eff Uc) 



:= C 2 o;o(^)- 



1 + 2mC 2 o;i(a;) 

Therefore, for low energy particles where one may neglect terms involving 4th and higher 
powers of p, the Hamiltonian h and consequently H describe motion of a particle with an 
effective position dependent mass m e g(x) that interacts with the potential w(x). Figure 5 
shows a graph of m e s (x) for m = l/2,h=l,L = 2 and ( = 1/3. For the same values of these 
parameters, w(x) = a (x)/9. See Figure 4 for a graph of a . 

If we replace (x,p) of (JHHj) and (fMj) with their classical counterparts (x c ,p c ), we obtain the 
'classical' Hamiltonian: 



2 a2 x 



2m 2 ^ nl c ^ c lS ; 2m eff (x c ) 



2 *2 00 

''" + w(x c ) + ^Y1 a ^ + °^ 3 )' ( 85 ) 

n=2 
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Figure 4: Graph of a , a x and a 2 
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Figure 5: Graph of the effective mass m e g (full curve) for m 
\, h = 1, L = 2 and ( = |. The dashed curve represents m = \ 



which coincides with the free particle Hamiltonian outside the physical interaction region, i.e., 
( — iT'ir)- The fact that this region is three times larger than the support (—-§,10 of the 
potential v(x) is quite surprising. Note also that H c is an even function of both the position x c 
and momentum p c variables. 

Figure 6 shows the phase space trajectories associated with the Hamiltonian H c for L = 
2, h = 1, m = 1/2, ( = Z = 1/3. For large values of the momentum the trajectories 
are open curves describing the scattering of a particle due to an interaction that takes place 
within the physical interaction region, (—3,3). For sufficiently small values of the momentum 
closed trajectories are generated. These describe a particle that is trapped inside the physical 
interaction region. This is consistent with the fact that for small p c , H c is dominated by the 
potential term w(x c ) which in view of its relation to ao(x) and Figure 4 can trap the particle. 

We wish to emphasize that because we have not yet takes the h — > limit of H c , we cannot 
identify it with the the true classical Hamiltonian H c for the quantum Hamiltonian h and 
consequently H. Given the limitations of our perturbative calculation of H c , we are unable to 
determine this limit. 9 Therefore, we cannot view the presence of closed phase space trajectories 
for H c as an evidence for the existence of bound states of h and H. This is especially because 
these trajectories are associated with very low momentum values where the quantum effects 
are expected to be dominant. 



5 Conclusion 

In this paper we explored for the first time the utility of the methods of pseudo- Hermit ian 
quantum mechanics in dealing with a non-Hermitian PT-symmetric potential v (x) that has 
a continuous spectrum. Using these methods we were able to obtain the explicit form of the 

9 This is in contrast with both the 'PT-symmetric square well and the T^T-symmetric cubic anharmonic 
oscillator studied in and , respectively. In the former system the presence of an exceptional spectral point 
imposes the condition that £ must be of order H 2 or higher and consequently the classical system is the same as 
that of the Hermitian infinite square well 5 . In the latter system, the h — + limit of the associated Hermitian 
Hamiltonian can be easily evaluated and classical Hamiltonian obtained 
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Figure 6: Phase space trajectories of the Hamiltonian H c (x c ,p c ) for 
m — 7},h — 1, L — 2 and C — |- The horizontal and vertical axes 
are respectively those of x c and p c . 

metric operator, the pseudo-Hermitian position and momentum operators, the localized states, 
and the equivalent Hermitian Hamiltonian perturbatively. 

Our analysis revealed the surprising fact that the physical interaction region for this model 
is three times larger than the support of the potential, i.e., there is a region of the configuration 
space in which v(x) vanishes but the interaction does not seize. 

A simple interpretation for this peculiar property is that the argument x of the potential 
v(x) is not a physical observable and the support (—■§,-§) of v(x) being a range of eigenvalues of 
x does not have a direct physical meaning. This observation underlines the importance of the 
Hermitian representation of non-Hermitian (inparticular PT-symmetric) Hamiltonians having 
a real spectrum. 

The Hermitian representation involves a nonlocal Hamiltonian that is not suitable for the 
computation of the energy spectrum or the ^-matrix of the theory. Yet it provides invalu- 
able insight in the physical meaning and potential applications of pseudo-Hermitian and VT- 
symmetric Hamiltonians and is indispensable for the determination of the other observables of 
the corresponding quantum systems. 
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